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Abstract 


Optimal control theory is employed to determine the performance of abort to 
orbit (ATO) and return to launch site (RTLS) maneuvers for a single-stage to orbit 
vehicle. The vehicle configuration examined is a seven engine, winged-body vehicle, 
that lifts-off vertically lands horizontally. The abort maneuvers occur as the vehicle 
ascends to orbit and axe initiated when the vehicle suffers an engine failure. The 
optimal control problems are numerically solved in discretized form via a non-linear 
programming (NLP) algorithm. A description highlighting the attributes of this NLP 
method is provided. 

ATO maneuver results show that the vehicle is capable of ascending to orbit with 
a single engine failure at lift-off. Two engine out ATO maneuvers are not possible 
from the launch pad, but are possible after launch when the thrust to weight ratio 
becomes sufficiently large. Results show that single engine out RTLS maneuvers can 
be made for up to 180 seconds after lift-off and that there are scenarios for which 
RTLS maneuvers should be performed instead of ATO maneuvers. 
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Chapter 1 


Introduction 


In 1993, Congress asked the National Aeronautics and Space Administration (NASA) 
to address the future needs of U.S. space transportation systems. Although current 
systems meet the functional needs of the U.S. space program, the systems lack cost 
effectiveness, reliability and operability. As a result, the ability of U.S. industry to 
compete in the international launch market has suffered. 

In response to the Congressional request, NASA studied [1] possible architectures 
for the next generation of launch vehicles. This next generation of vehicles is to ad- 
dress the deficiencies of the current systems and at the same time focus on increasing 
safety. Recognizing the major advances in technology and vehicle design that have 
been made within the past decade, the NASA recommended the development of a 
fully reusable single-stage to orbit (SSTO) vehicle. 

Addressing the issue of increased safety, the purpose of this work is to determine 
the abort performance of abort to orbit (ATO) and return to launch site (RTLS) 
maneuvers for a SSTO vehicle. The vehicle configuration considered is a seven en- 
gine, winged-body vehicle that lifts-off vertically and lands horizontally. The abort 
maneuvers occur as the vehicle ascends to orbit and are initiated when the SSTO ve- 
hicle suffers an engine failure. A point mass model is used to describe the motion of 
the SSTO vehicle. A non-linear programming (NLP) method is developed for solving 
optimal control problems is used in the analysis. 

Before the abort maneuvers are studied, a baseline ascent trajectory is generated. 



In this ascent the vehicle has a full complement of engines. The trajectory is shaped 
so as to maximize the amount of mass delivered to orbit. This trajectory provides a 
reference to which the abort maneuvers can be compared. 

The amount of mass delivered to orbit is also used as a performance measure for 
the ATO maneuvers. The performance for problems with various engine failure times 
is obtained. The vehicle travels on the baseline ascent trajectory until the time of 
engine failure. Single-engine and two-engine-out scenarios are examined. 

Several problems are investigated for the RTLS maneuvers. Like the ATO ma- 
neuvers, the vehicle initially travels on the baseline ascent trajectory before suffering 
an engine failure. The maximum time after launch that the vehicle can perform an 
RTLS maneuver is determined. Also, the limits on the payload mass that can be 
returned to the launch site are determined for specified engine failure times. 

Previously [2], abort capabilities for a similar vehicle were determined. In this 
study, the maximum allowable normal force and terminal altitude for the RTLS ma- 
neuvers are significantly larger than those allowed for this study. 

The results of computations for a winged-body SSTO configuration are given in 
this paper. Also, the attributes of the NLP method are described. 
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Chapter 2 


SSTO Background 


A single-stage launch vehicle discards only propellants while ascending to orbit. For 
a single-stage vehicle, each pound of structure trades off for a pound of payload; 
therefore minimizing the dry weight of the vehicle is crucial in the design process. 
Because technologies required to produce adequate performance margins have been 
lacking, the feasiblity of a single-stage launch vehicle has been questioned. 

In constrast, a multiple-stage launch vehicle discards both propellants and struc- 
tural weight while ascending to orbit. By staging, the weight and performance re- 
quired by each stage is reduced, thus making the multiple-stage launch vehicle a 
feasible alternative. However, staging also introduces undesirable operational com- 
plexities and cost. 

Recent technological advances in the fields of material science and propulsion have 
influenced NASA to propose the development of a single-stage launch vehicle. One 
of the proposed designs is a winged-body configuration. 

The selection and design of this configuration is described in [3]. As illustrated 
in Figure 2.1, this vehicle is powered by seven space shuttle main engine (SSME) 
derivative engines. This vehicle launches vertically, and lands horizontally. The 
vehicle uses chemical propulsion with liquid hydrogen (LH) and liquid oxygen (LOX) 
serving as the fuel and oxidizer respectively. As shown in Figure 2.1, the payload bay 
separates the forward LOX and aft LH tanks. In table 2.1 various vehicle parameters 
are given including the performance characteristics of the SSME derivative engines. 
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Wing Reference Area = 4192 ft 2 



Payload Bay : 

Diameter = 15 ft, Length = 30 ft 

Figure 2.1: Winged-Body Configuration 

Issues that must be addressed in order for the SSTO vehicle to demonstrate im- 
proved operability are treated in [4]. This document discusses lessons that have been 
learned from the space shuttle program an other launch vehicle concepts. 

In [5] requirements for the SSTO vehicle are outlined. This study addresses the 
requirement that the vehicle is to have “At minimum, single-engine-out capability 
throughout ascent (either return to site or abort to orbit/once around)”. 
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Table 2.1: Miscellaneous Vehicle Parameters 


Number of Engines 

7 

Vacuum Thrust 

463.90 

klbf 

Vacuum Specific Impulse 

447.3 

sec 

Engine Mass 

6.79 

klbm 

Oxider To Fuel Ratio 

6.0 

Exit Area 

27.55 

ft 2 

Vehicle Dry Mass 

206.50 

klbm 

Vehicle Gross Mass at Lift-Off 

2383.43 

klbm 

Nominal Payload Mass 

20.0 

klbm 

Vehicle Reference Area 

4192.2 

ft 2 


5 






















Chapter 3 


SSTO Model Description 


In this section, the physical model used for the abort maneuvers is discussed. Equa- 
tions describing the motion of the vehicle, the acceleration due to gravity, the atmo- 
sphere, the aerodynamics and the propulsion are given. 


3.1 Equations of Motion 

The center of mass motion for a thrusting SSTO vehicle flying over a rotating Earth 
can be described by (see appendix A for development) 

. > 

6 = v cos cos rp / (r cos p) 

<p = ucos 7 sin ip/r 
h = usin7 

v — (Tcosa — D)/m — g sin 7 + u 2 r cos <p(sm 7 cos <j> — cos 7 sin <f> sin ip) 
ip = (T sin a -|- L) sin g/(mv cos 7) —u> 2 r sin <p cos <p cos tp/(v cos 7) ' 

+2u;(tan 7 cos <j> sin ip — sin^)/(ucos7) — ( v/r ) cos 7 cos ip tan <f> 

7 = (Tsino -f L) cos fi/(mv) + (u> 2 r/v) cos <?i(cos 7 cos ^ + sin 7 sin^sin^) 

-f 2 uj cos <j> cos ip + ( v/r — g/v)cos~f 

ril = N en gines^lT va cJ sp) J 

( 3 . 1 ) 

The state variables for this system of equations are longitude 8, latitude <j>, altitude 
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h, airspeed u, relative heading 0, relative flight-path 7 and mass m. The controls are 
attack angle a, bank angle fi and throttle r}. These equations assume that moment 
equilibrium can be maintained by the control system. These equations also assume 
that a positive bank angle generates a heading to the north for a vehicle flying west 
to east. 


3.2 Gravity Model 

The Earth is assumed to be a sphere whose radius r e represents mean sea level. If 
r is the radial distance from the center of the Earth to the center of gravity of the 
vehicle, the acceleration of gravity is given by the inverse square law 

9 = go{r E /r ) 2 (3.2) 


3.3 Atmosphere Model 

A simple exponential model was adopted for the atmosphere model. In this model 
the atmospheric pressure is calculated by 


P = Poe 


-WPp) 


Similarly, the atmospheric density is found with 


(3.3) 


P = pot 


-(»■/».) 


(3.4) 


Assuming that air behaves as an ideal gas, the temperature is calculated using 


T = p/ipR,,,,) 


(3.5) 


Table 3.1: Atmospheric Parameters 


TE 

2.0925721e+07 

ft 

Po 

2.1162e+03 

psf 

Pv 

2.3127e+04 

ft 

Po 

2.3769e-03 

slugs/ft 3 

Pr 

2.3341e+04 

ft 

Rgas 

1.7163e+03 

ft 2 /sec 2 °R 

k 

1.4 


Since air is assumed to behave as an ideal gas, the speed of sound may be calculated 
with 



(3.6) 


Table 3.1 summarizes the values of the constants used to evaluate atmospheric quan- 
tities. 


3.4 Aerodynamics 

The lift and drag acting on the SSTO vehicle are related to the lift and drag coefficients 
(cl and cp) as follows 


L = ~pv 2 Sc L 

(3.7) 

D = ^ pv 2 Sc D 

(3.8) 


where S is the aerodynamic reference area. The lift and drag coefficients are functions 
of mach number and angle of attack. Mach number is found with 

M =- (3.9) 

a 
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The aerodynamic coefficients are modeled with 


c L (M,a) = c Lq (M) a 


(3.10) 


c d (M, a) = c Do (M) + c Da (M) c* 2 (3.11) 

In this model, the coefficients ci a and Cu Q are chosen so as to minimize the Euclidean 
norm of the error between the model and existing tabular data for the lift and drag 
coefficients at the table values. The Euclidean norm of the error between the lift 
model and the lift coefficient data ( cltab ) is described by 


= \i C La Cij - c LTAB {aj )] 2 - - < a, < ^ (3.12) 

Note that only the tabular lift coefficients corresponding to attack angles in the 
interval [ — 7r/2, tt/ 2] influence Jl- The value of ct Q that minimizes Jl is 


a j c Ltab ( a i ) 

L ° = EX 


7T 7 r 

— — < a, < — 
2 ~ 1 ~ 2 


(3.13) 


Similarly, the Euclidean norm of the error between the drag model and the drag 
coefficient data ( cdtab ) is described by function 


7T 7T 

Jv = Y! 2 ^ CDo + CDa “ C DTAB(<*j )] 2 ~ ^ < Oij < - 

3 


(3.14) 


where, for simplicity, the coefficient c Bo was chosen to be the value of the tabular 
drag coefficient for zero angle of attack. The value of c Ba that minimizes Jv is 


T,j[cDTABi a j) ~ c DoWj 

CDa = 

^3 3 


7T 7T 

— < a, < — 
2 - 3 ~ 2 


(3.15) 


Figures 3.1 through 3.4 illustrate the performance of the aerodynamic models in 
comparison to the tabular data. Surface plots of the aerodynamic models for the lift 
and drag coefficients are provided in Figures 3.1 and 3.2. In these figures the tabular 
data points are included as asterisks to provide an idea of how accurately the models 
approximates the tabular data. The shapes of these surfaces in the transonic region 
make modeling the surfaces with a simple analytic expression difficult. 
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Tabular Lift Coefficients 
Modeled Lift Coefficients 





Defining the error in the aerodynamic models as the difference between the mod- 
eled coefficients and the tabular coefficients, Figures 3.3 and 3.4 illustrate the absolute 
error in the aerodynamic models for the lift and drag coefficients respectively. The 
tabular data consists of three different subtables corresponding to subsonic, super- 
sonic, and hypersonic mach ranges. Each of the subtables are defined with different 
attack angles. Because the tabular data consists of three subtables, three different 
surfaces are required to describe the absolute error. These figures show that the 
aerodynamic models adequately approximate the tabular data. 



Figure 3.3: Comparison of Model and Tabular Lift Coefficients 
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Figure 3.4: Comparison of Model and Tabular Drag Coefficients 


Values for CL a , cd 0 and cd 0 are defined for various mach numbers in table 3.2. 
Interpolation is used to determine the values of these coefficients for mach numbers 
not represented in the table. The interpolation technique used to evaluate these 
coefficients is an alternative method to piecewise linear interpolation. 

Many optimization schemes, including the one employed in this study, require 
derivative information. In general, piecewise linear interpolation does not provide 
first derivative continuity. For this reason, optimization schemes that use piecewise 
linear interpolation may fail. The technique used in this analysis guarantees second 
derivative continuity. A detailed description of this interpolation technique is given 
in [6] 
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Table 3.2: Coefficients Used in Aerodynamic Models 


Mach 

P Lot 

CDo 

P Da 

0.0000 

2.5004 

0.0277 

1.1678 

0.3000 

2.5004 

0.0277 

1.1678 

0.6000 

2.6467 

0.0388 

1.2819 

0.9500 

2.0782 

0.0204 

1.3840 

1.1000 

2.4824 

0.4071 

3.9257 

2.0000 

1.6739 

0.0918 

2.6077 

3.0000 

1.3597 

0.0811 

2.1939 

4.0000 

1.2436 

0.0767 

2.0492 

6.0000 

1.1742 

0.0727 

1.9773 

8.0000 

1.1465 

0.0710 

1.9499 

10.0000 

1.1335 

0.0709 

1.9389 

12.0000 

1.1497 

0.0720 

1.8792 

15.0000 

1.1402 

0.0704 

1.8711 

20.0000 

1.1285 

0.0710 

1.8550 

25.0000 

1.1208 

0.0739 

1.8314 


3.5 Propulsion 

This SSTO vehicle model has seven Space Shuttle Main Engine (SSME) derivative 
engines. The total thrust of the SSTO vehicle is given by 


T — A engines {V^vac A e p) 
and pertinent engine parameters are found in table 2.1. 


(3.16) 
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Chapter 4 


Solution Method 


This section begins by discussing optimal control at a generic level. Single phase and 
multiple phase optimal control problems are described. Next, a numerical approach 
for transforming the optimal control problems into algebraic non-linear programming 
(NLP) problems is outlined. Finally, the performance of the algorithm is illustrated 
by solving a simple example problem. 


4.1 Optimal Control 

The objective of optimal control theory is to determine the control signals that will 
cause a system to satisfy the physical constraints and at the same time minimize some 
performance criterion. 

4.1.1 Single Phase Optimal Control Problem 

For a specified initial time f 0 , choose the initial states :r(fo) € 7l n , final time f/, 
control history [u(t) £ 'R m t 0 <r <tf] and free parameters z € Wto minimize the 
scalar objective function 


\J ^[^(^0)5 


(4.1) 
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The objective function or performance index is subject to the dynamical con- 
straints 

x (t) - (x(t 0 ) + /[x(r),u(r), 2 ]dr^ = 0 (4.2) 

and boundary conditions specified as 

W[x(to),x(t f ),t f ,z] = 0 (4.3) 

The control function u(t) is assumed to be integrable on the interval t 0 < t <tf. 

Inequality constraints may also exist. These constraints appear in either of two 
forms: 

state/control 


or state-only 


C(x, z) < 0 


(4.4) 


S(x,t,z)< 0 (4.5) 

The constraints are categorized in this manner because the control of S(x,t,z) is 
obtained only by changing the qth time derivative of the constraint S q (x,t,z) which 
depends explicitly on the controls u. As the state trajectory enters the constraint 
boundary <S(x,f,z) = 0, the tangency constraints 

[ S(x,t,z ) = 0 
x,t,z) = 0 

[ <S ?-1 (x,f,z) = 0 

must be satisfied. These tangency constraints apply for state trajectories leaving the 
constraint boundary as well. Since the control appears explicitly in 4.4, “tangency” 
is not required for these inequality constraints [7]. As a result satisfying state/control 



(4.6) 
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inequality constraints is typically somewhat easier than satisfying state-only inequal- 
ity constraints. 

A useful example of a free parameter 2 is the duration of the trajectory in a free 
time problem. The state equations for a free time problem 


— f(x,u,t) t 0 <t<tf t 0 fixed and t/free 
at 


can be recast as 


— = {t s - t 0 )f(x, u,t) 0 < r < 1 


where 


r-(t-to)/(tf-t 0 ) (4.9) 

One technique for solving a free time problem is to select a free parameter to represent 
the duration of the trajectory tj — to. An additional inequality constraint 


tj - t 0 > 0 


(4.10) 


prevents negative time scaling of the state equations thus ensuring that time mono- 
tonically increases. 

4.1.2 Multiple Phase Optimal-Control Problem 

Single phase optimal control problems may only have boundary conditions at the 
initial and final times of the trajectory. Frequently, a need arises for additional 
boundary conditions in the interior of the interval [t 0 ,tf\. 

Multiple phase optimal control problems allow for the existence of interior bound- 
ary conditions. The states and controls of each phase of a multiple phase problem 
can be represented by 


Xk k = 1, . . . ,L Xk € TZ nk 
Uk k = l,...,L Uk£'R. mk 


(4.11) 
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where L is the total number of phases. If the initial and final times of each phase are 
denoted by tko and tkj respectively, the initial and final values of the states of each 
phase are 

Xko = %k(tko') k 1, . . . , X 
Xkf = %k(tkf') k = 1* • • • j L 
Figure 4.1 shows a time-line for a multiple phase problem. In this figure, the phases 
are denoted by 


( 4 . 12 ) 


*1 *2 


* * 


b 


t 


t 


+ 


Figure 4.1: Multiple Phase Time-line 
The objective function for a multiple phase problem is of the form 

i 7 = J[xiO, Xif, X20i X2 ft ' ' * 5 Xjj 0, X£,y, 2] 

The multiple phase dynamical constraints are 


(4.13) 


x k (t) - (x k (tk 0 ) + J t f[xk(r), Ufc(r), z\drj =0 k = l,...,L (4.14) 


with the boundary conditions 


\H[xio, Xi y, X20; X 2 y, ■ • • , Xlq , Xlj , z] — 0 


( 4 . 15 ) 
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The inequality constraints are expressed as 


C(xk,Uk,t,z) < 0 k = 1,...,L (4.16) 

and 

S(xk,t,z) < 0 k = l,...,L (4-17) 

4.2 Numerical Approach 

In this approach, optimal control problems are transformed into algebraic NLP prob- 
lems of the form 


minimize J{y) 
yen M 

subject to : 


(4.18) 


E(y) = 0 

i{y) > 0 ^ 

where y is a set of design variables used to define a scalar objective <f>{y] that is subject 
to a set of equality constraints E(y) and a set of inequality constraints I(y). The 
design variables y will be defined shortly. 

For the sake of simplicity, let us first consider a single phase optimal control prob- 
lem with free parameters. This approach assumes that the duration of the trajectory 
is unity. Often, this requires a scaling of the state equations using the same method 
previously illustrated by equations 4.7 through 4.9. 

The transformation process begins with the discretization of the state equations 
into N equally spaced intervals resulting in N+l nodes. Since the duration of the 
trajectory is unity, the time step between nodes is 



(4.19) 
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Therefore, the time at each of the nodes is given by 


= + 1 (4-20) 

The design variables y to be optimized in this approach are the values of the 
states X{ € TIT at the N+l nodes, the values of the controls u t - 6 7i m at the N nodal 
midpoints and the free parameters z € TV . The ordering of the design variables is 

y = [xi u A x 2 u 2 x N u N x N+1 z] T (4.21) 

The objective function to be minimized for a single phase problem is given by 

J = J{xi,x N+ i,z) (4.22) 

Using a second-order, mid-point integration rule, the state equations produce the 
following system of equality constraints 

x i+i f (®» j Wj, z'j — 0, 2 = 1 ,..., N (4.23) 

where 

S,- = |(x, + x i+ i) (4.24) 

Another set of equality constraints are the boundary conditions for the state equations 
which are of the form 


®(x l5 a;jv + i,2) = 0 
State/control inequality constraints are of the form 


C,(xi,ui,z) 

Ci(x 2 jU 2: z) 


> 0 i = 1, .. .,n c 


Ci(x N , u N , z) 


(4.25) 


(4.26) 
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where nc is the number of state/control constraints. Like the discretized state equa- 
tions, the state/control inequality constraints are enforced at the N nodal mid-points. 
The state-only inequality constraints are of the form 

<S t (xi,z) 

Si(x 2 , z) 

Si (xjv+1 , zj 

where n, $ is the number of state only constraints. Finally, free parameter inequality 
constraints expressed as 


>0 * = 1,..., (4.27) 


V{z) > 0 (4.28) 

may also exist. 

When integral cost expressions such as 

J = [ C(x, u, z)dr (4.29) 

J o 

must be addressed, the cost function is defined as a state. The state equation for this 
new state is 


xj = C(x,u,z ) (4.30) 

Now consider a multiple phase optimal control problem. The state equations for 
an L phase problem are expressed as 

Ft = fi :(xfc, Ufc, x) k — 1, . . . , L (4.31) 
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where Xk € TZ nk , Uk 6 'K mk and z E W . After first scaling these equations such that 
the duration of each phase is unity, the equations are discretized into Nk equally 
spaced intervals. The NLP problem is formulated as in 4.18, except the design vari- 
ables y are now ordered as 

y = [y\ vi Vl zf (4.32) 

where the design variables for the kth phase are 

Vk [ x lk k x 2k ^2 k %Nk k 'U'Nkk (4.33) 

Denoting the initial and final values of the states for the kth phase by xi* and x^ kk 
respectively, the performance index is expressed as 

3 kJ [®11j x Nhi x N22’> ' ‘ * j x \ Li z] (4.34) 

Using the same mid-point integration rule, the state equations yield the following set 
of equality constraints 

x t+11 - x tl - f(xn, uh,z)/N i = 0; * = 1, . . . , JVi 
x, + i 2 - x,- 2 - f(xi 2 , u i2 , z)/N 2 = 0; i — 1,..., N 2 

Xi+lL - x iL - f&L, U iL , z)/Nl = 0 ; i = 1 , . . . , N l 
The boundary conditions now take the form 

^[*11, x Nni ^12? X N^2 , • ’ ’ j XlL, Xn l i,i z] — 0 (4.36) 

The state/control inequality constraints for the kth phase are described as 

Uj(xi ki Ui ki z) 

Ci(x 2k ,U 2 k,z) 

Ci(x Nkk ,u Nkk ,z) 
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Similarly, the state-only inequality constraints for the kth phase are 


Si(xi k ,z) 

Si(x 2k ,z) 


>0 i = l,...,ncfc 


(4.38) 


[ Si(x Nk+lk ,z) J 

The inequality constraints are concatenated in the same manner as the state equation 
equality constraints yielding the two sets of inequality constraints 



"Cl ' 


' Si ' 

c = 

c 2 

S = 

S 2 


Cl _ 


S L 


(4.39) 


Free parameter inequality constraints of the form 4.28 may also be included. This 
approach is not self-starting. An initial guess is required to start an iterative solution 
process. The Fortran code NPSOL [8] was selected to solve the NLP problem. In 
appendix B, the codes that provided the skeleton of this method are described in 
detail. 


4.3 Algorithm Performance 

As a brief illustration of the algorithmic performance of NPSOL and the discretized 
problem representation described above, a simple problem will be solved. The solu- 
tion will then be compared to the closed-form solution. 

The double integrator problem is posed as follows : 

Minimize the cost function 



(4.40) 
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subject to the differential constraints 


i\ = x 2 ; £i(0) = 1 *€[0,1] 

x 2 = u; x 2 (0) = 0 

and the boundary conditions or terminal constraints 


(4.41) 


4>i : Xi(l)-1 = 0 

► 

$2 : x 2 (l) = 1 

The integral cost function is addressed as in 4.29 and 4.30. 
is 


(4.42) 


The resulting cost function 


J = ® 3 ( 1 ) 

which is subject to the state equations 


(4.43) 



Xl 

= x 2 ] 

xi(0) 


X2 

= u; 

x 2 (0) 


x 3 

= lu 2 

xs(0) 

the terminal constraints 

remain the 

same. 

is 






Xl (<) = 

3 1 3 - 



x 2 (t) = 

9 1 2 - 



Xs(t) = 

54 * 3 



u(i) = 

18* - 


(4.44) 


(4.45) 


The optimal value of the objective function is 14. 
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Groundtrack Comparison 



Figure 4.2: Ground-track Comparison 

Solutions for this problem were also found for discretization intervals of N = 
(5,10, and 20). Figure 4.2 shows the ground-tracks of the numerical solutions in 
comparison to 4.45. In table 4.1 results from the three different nodal densities studied 
are compared. As the example problem illustrates, the accuracy of this algorithm 
is enhanced by increasing the nodal density. As expected, as the nodal density is 
increased, the computational effort also increases. Thus, a trade-off must be made 
between desired accuracy and exceptable computational effort. 


Table 4.1: Numerical Approach Results 


Discretization Intervals 

5 

10 

20 

50 

Design Variables 

23 

43 

83 

203 

Iterations 

15 

24 

34 

27 

Cost Function Evaluations 

20 

55 

44 

39 

CPU Time, sec 

0.5 

2.7 

20.4 

203.1 

Optimal Cost 

14.5625 
■ 

14.1364 

14.0338 

14.0054 
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Chapter 5 


Abort Maneuvers 


In this section the abort capabilities of the SSTO vehicle are examined. Each of the 
optimal control problems used in this investigation are described. Results for each 
problem are presented after the problem is described. 


5.1 Baseline Performance Problem Description 

Before examining the abort maneuvers, a baseline for performance must be created. 
For this study, the baseline will be a maximum mass to orbit ascent trajectory for a 
vehicle with a full complement of engines. The target orbit is a 50 by 100 nmi orbit 
with an inclination of 51.6°. 

If m/ is the mass of the vehicle at orbit insertion, the baseline trajectory is found 
by minimizing the objective function 


3 = -m } (5.1) 

Note that minimizing — mj is equivalent to maximizing m/. The state equations 
for this problem were previously described in 3.1. The boundary conditions for this 
problem are 
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Vo 

— 

0.1 

ft/s 


ho 

= 

0 

ft 


m 0 

= 

2383.43 

klbm 


0o 

= 

-80.71 

deg 


<f>0 

= 

28.47 

deg 

> 

Vf 

= 

25852 

ft/s 


7/ 

= 

0 

deg 


hf 

= 

3.03805 

ft 


if 

= 

51.6 

deg 


As shown in [11], the final inclination if can be found with 


(5.2) 


VfCOS~ifCos(f>jCos ij)j + u>(hj + r e )cos 2 <t> 
sjiyfcos^jcos'ipj + ( hf + r e )cocos(f>) 2 + (yfCos'ffsim^f) 2 

The initial flight-path and heading angles are designated as free variables, thus ne- 
glecting the brief and highly constrained motion required to clear the launch pad. 
The state/control inequality constraints for this problem are 




0.25 < T) < 1.09 ' 

\a\ < 45° 

(5.4) 

\F n \ < 3.81232e + 05 lbf 
\A\ < 3<7o 

where A is the acceleration of the vehicle, g 0 denotes the gravitational acceleration 
at sea level, and Ffj is the normal force, given by 


F n = 1/2 pv 2 Sc N 

where the normal force coefficient cjv is 


(5.5) 


cjv = cd sin oc + cjj cos a 


(5.6) 
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The state-only inequality constraints are 


h 

> 

0 

ft 

m 

> 

TTldry 


q 

< 

1000 

psf 

M 

> 

0 


ItI 

< 

90° 



where the dynamic pressure q is given by 


(5.7) 


q = l/2pv 2 (5.8) 

Since this is a free time problem, the duration of the trajectory r is a free param- 
eter. The state equations 3.1 must be scaled by this parameter. The free parameter 
inequality constraint shown below guarantees that time will monotonically increase. 


r > 0 


(5.9) 


5.2 Baseline Performance Results 

The resulting state trajectories for the baseline performance problem are shown in 
Figures 5.1 through 5.6. The mass lifted into orbit is 286.74 klbm which is 12 percent 
of the mass of the vehicle at lift-off. If one assumes that all of the propellant is used 
during the ascent, the corresponding maximum payload is 80.25 klbm. The time of 
flight is 370 seconds. 
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Figure 5.1: Baseline Ascent Trajectory : Air Speed 



Figure 5.2: Baseline Ascent Trajectory : Relative Flight-Path 
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Altitude, ft 



Figure 5.3: Baseline Ascent Trajectory : Altitude 



Figure 5.4: Baseline Ascent Trajectory : Mass 


29 




3i 




The control histories for this ascent are illustrated by Figures 5.7 through 5.9. 
The throttle history shows that the vehicle flies full throttle for most of the first 
160 seconds of the ascent. At this time the vehicle begins to reduce the throttle 
to prevent the violation of the acceleration limit. A slight reduction in throttle is 
seen 50 seconds after launch. This reduction is to prevent the violation of the normal 
force constraint. Figure 5.10 shows when various inequality constraints become active 
during the ascent. In this figure, active constraints are represented with lines passing 
through asterisks. 



Figure 5.7: Baseline Ascent Trajectory : Throttle 
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Figure 5.8: Baseline Ascent Trajectory : Attack Angle 



Figure 5.9: Baseline Ascent Trajectory : Bank Angle 
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Active Inequality Constraints vs. Time 


3 
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0 50 100 150 200 250 300 350 400 

Time, sec 

Figure 5.10: Active Inequality Constraints in Baseline Ascent 

In Figure 5.11 the attack angle is plotted as a function of mach number. This figure 
shows that as the vehicle travels through the transonic mach region the magnitude 
of the attack angle is reduced. Figure 5.12 shows the drag profile for the trajectory. 
The maximum value of drag is 845 klbf and occurs 70 seconds after lift-off when 
the vehicle is traveling at mach 1.18. Since the drag acting on the vehicle reaches a 
maximum in the transonic region (denoted by asterisks), the vehicle must fly with 
low angles of attack in order to minimize drag and increase performance. 
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Attack Angle vs. Mach 







5.3 ATO Maneuver Problem Descriptions 


In this problem, the vehicle ascends to orbit by traveling along the previously de- 
scribed baseline trajectory. At time after launch, the vehicle suffers an engine 
failure. The objective of this problem is to determine the loss in performance caused 
by an engine failure at time t^. Single-engine and two-engine-out scenarios are 
studied. 

As in the baseline trajectory problem, the performance index for this problem is 

J = —mj (5.10) 



The state equations 3.1 are modified to account for the appropriate reduction in 
the number of engines N eng i nes . In Figure 5.3 provides an illustration of the ATO 
maneuver problem. In this figure, the baseline state trajectory is denoted as £s(f). 
Refering to this figure the boundary conditions for this problem are 
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Uo 

= 

vb^a) 


7o 

— 

1b{Ia) 


ho 

= 



m 0 

= 

m B {t A ) 


*l > 0 

= 



00 

= 

0b(Li) 


4>o 

— 

^b^a) 


v f 

= 

25852 

ft/s 

7/ 

= 

0 

deg 

hf 

— 

3.03805e + 05 

ft 

*/ 

= 

51.6 

deg 


(5.11) 


Since the baseline ascent trajectory is represented by a discrete number of states, 
interpolation is used to determine the values of xb^a) for which data points do not 
exist. 

The state/control and state-only inequality constraints are the same as those 
enforced in the baseline trajectory previously given in 5.4 and 5.7. Like the the 
baseline trajectory problem, this problem is a free time problem. As always, the state 
equations must be scaled by the duration of the trajectory r which is represented by 
a free parameter. The free parameter inequality constraint 


r > 0 


(5.12) 


is added to prevent the time scale from become negative. 


5.4 ATO Maneuver Results 

Figure 5.13 shows the maximum payload mass that can be delivered to orbit as a 
function of the time of engine failure t^. Results from single-engine and two-engine- 
out scenarios are shown. 
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In this figure, initial results for the single engine out scenario are labeled Single- 
Engine-Out 1. These results were obtained as was varied from 0 to 100 seconds 
after launch in 10 second intervals. As t a was varied, the solution for the problem 
with current value of t a served as the initial guess for the next problem with 
incremented. These results suggest that a slight decrease in performance occurs as 
tA is increased from 50 to 60 seconds after launch. The large increase in performance 
shown as t.A is varied from 80 to 90 seconds after launch hints that the “solutions” 
for engine failure times in the interval 60 < < 80 are local minima. 

The Single-Engine-Out 2 results were obtained by walking the time of engine 
failure from 90 to 40 seconds after launch in 5 second intervals. Again, each solution 
served as the initial guess for the problem with the subsequent value of The 
improved performance shown by the Single-Engine-Out 2 solutions supports the 
previous intuition that the Single-Engine-Out 1 solutions for this same time range 
are local minima. Note that since Single-Engine-Out 2 results were not determined 
to be the global minimum (a difficult task for most optimal control problems), these 
results should also be categorized as local minima. Both sets of solutions are feasible 
solutions. 

Solutions for two-engine-out cases were found for engine failure times in the in- 
terval 50 < tA < 100. Solutions could not be found for engine failure times below 50 
seconds after launch. If traveling on the baseline ascent trajectory, a vehicle with two 
failed engines has a thrust to weight ratio less than 1 until 25 seconds after lift-off. 
The thrust to weight ratio remains inadequate for the two-engine-out ATO maneuvers 
until the conditions 50 seconds after launch. 
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Maximum Mass of Payload vs. Time of Engine Failure 



Further investigation of the single-engine-out results shows that a decrease in 
performance occurs as the time of engine failure varies from 50 to 60 seconds after 
launch. As Figure 5.12 previously illustrated, the vehicle begins to enter the transonic 
region during this time frame. Figure 5.14 shows the two drag profiles for these to 
single-engine-out ATO maneuvers. This figure shows that the 60 second case flies 
with more drag than the 50 second case resulting in the decreased performance. 
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Drag vs. Time : Single Engine Out 



Time, sec 

Figure 5.14: Drag Profiles for ATO Maneuvers Initiated at 50 and 60 Seconds After 
Launch 

Figure 5.15 shows that single-engine and two-engine-out ATO maneuvers initiated 
at 60 seconds after launch fly more lofting trajectories than the baseline trajectory. 
Note that the two-engine-out scenario flies a more lofting trajectory than the single- 
engine-out scenario. To investigate this behavior, examine the expression for 
approximated by 

Jl 

— = mv sin 7/(T cos a — D) (5.13) 

dv 

Because of the engine failure(s), T cos a is larger for the baseline ascent than for the 
ATO trajectories. As a result ^ is larger for the ATO trajectories. 
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Figure 5.15: Lofting of ATO Maneuvers Initiated 60 Seconds After Launch 

Figure 5.16 shows the attack angle profiles for the for the same two scenarios. This 
figure shows that the vehicle must fly with higher attack angles after engine failure 
occurs. 

Figure 5.17 illustrates mass as a function of time for the same scenarios. As 
expected the single-engine-out scenario out performs the two-engine-out scenario. 
Unsurprisingly, the time of flight for the two-engine-out scenario is longer than the 
single-out-scenario, since the vehicle has even less thrusting capability. 
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Attack Angle vs. Time : Abort Time = 60 



Figure 5.16: Attack Angle Profiles for ATO Maneuvers Initiated 60 Seconds After 
Launch 


Mass vs. Time 



Figure 5.17: Mass Profiles for ATO Maneuvers Initiated 60 Seconds After Launch 
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Finally, the throttle profiles depicted in Figure 5.18 illustrate that the reduction 
thrust due to engine failure is compensated for by maintaining the throttle at the 
maximum value longer. 


Throttle vs. Time 



Figure 5.18: Throttle Profiles for ATO Maneuvers Initiated 60 Seconds After Launch 
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5.5 RTLS Maneuver Problem Descriptions 


Several problems are investigated for the return to launch site maneuver. Each for- 
mulation consists of two phases and thus utilizes the multiple phase capability of 
the numerical approach previously described. Like the ATO maneuver problem, the 
RTLS maneuvers are initiated from fixed points on the baseline ascent trajectory. 
The first phase is designed to empty the propellant from the vehicle tanks. During 
the second phase, the lightened vehicle flies back to the landing site. Figure 5.19 
illustrates the different phases of the RTLS maneuver. 



Figure 5.19: RTLS Phases 
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5.5.1 Maximize Time of Abort 


The objective of this problem is to determine the maximum time after launch that a 
vehicle suffering from an engine failure can return to the launch site. This “point of 
no return” will be found for a vehicle carrying the nominal payload mass of 20 klbm. 
If the engine failure occurs at time t A after launch, the performance index for this 
problem is 


J = -t A (5.14) 

The state equations 3.1 are modified to reflect the reduction in the number of engines. 
If the baseline trajectory is denoted by x B (t ), the boundary conditions for this problem 
are 
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(5.15) 
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where the final conditions represent typical landing conditions. The state/control 
inequality constraints for this problem are 



0.25 

<t] < 1.09 

$// : 

0.0 

< ri < 0.0 ( i.e . rj = ( 


M 

^ CXmax 


\Fn\ 

< 3.81232e + 05 lbf 


\A\ 

< 3<?o 


(5.16) 


where the throttle limits are defined for each phase $ and the maximum attack angle 
tt mai £ [60°, 90°]. Note that in the second phase, the throttle is constrained to zero. 
Therefore, the thrust must also be set to zero or the back pressure term A e p will be 
incorrectly applied which may result in negative thrust. The state-only inequality 
constraints are 


h 

> 

0 

ft 

m 

> 

^ dry "f" YTlpayload. 


<? 
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1000 

psf 

M 
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Ob 
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< 


(5.17) 


where the mach number is constrained to be above the stall mach. 

Three free parameters exist for this problem: the abort time t a, the empty tank 
time ts and the terminal or landing time Denoting the final time of the baseline 
ascent trajectory as is/, the free parameter inequality constraints are 


0 < t A < tsf 
tE > l A 
tL > 


(5.18) 


5.5.2 Payload Returned To Launch Site 

The results from the ATO maneuver study have shown that the subject SSTO vehicle 
has single-engine-out capability from the launch pad. The trade-off is a reduction in 
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the amount of mass that can be delivered to orbit. Operationally, this dictates the 
scenario for performing an RTLS manevuer. A RTLS maneuver should be performed 
if the payload mass is greater than that which can be delivered to orbit by an ATO 
maneuver. This assumes that an RTLS maneuver can be made with such a payload 
mass. 

To test this assumption, the maximum and minimum amounts of payload that 
can be returned to the launch site are determined for specified times of engine failure 
tA ■ For the maximum payload problem, the objective function is 

\J — payload (5.19) 

while for the minimum payload problem 


«. y — ^payload (t).20) 

Again, the state equations 3.1 are modified to reflect the reduction in the number 
of engines. The boundary conditions, state/control inequality constraints and state- 
only inequality constraints are the same as those previously given by 5.15, 5.16 and 
5.17 respectively. 

In these problems the free parameters are the payload mass m pa yioad , the empty 
tank time t E and the terminal time t E . Recall, the baseline ascent trajectory delivers 
a maximum payload of 80.25 klbm to orbit. If the maximum payload is denoted as 
payload max , the free parameter inequality constraints for these problems are 


0 ^ m payload « P^ytoCt(I rrlax 

ts > tA 
tL > t E 


(5.21) 
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5.6 RTLS Maneuver Results 


Investigating a single-engine-out scenario, the maximum time of abort was determined 
for a vehicle carrying the nominal payload of 20 klbm. Various attack angle limits 
a max were examined. Figure 5.20 shows that the maximum time of abort decreases 
as the attack angle is more tightly constrained. 



Figure 5.20: RTLS Maximum Abort Time Versus Maximum Attack Angle For Single 
Engine Out Scenario 

To make the turn back toward the landing site, the vehicle must travel at reduced 
speeds to prevent exceeding critical structural loadings. The speed of the vehicle 
can be slowed by either decreasing the throttle or by increasing the drag acting 
on the vehicle. A vehicle flying with minimum throttle can only reduce speed by 
increasing drag which calls for high attack angles. As the attack angle is more tightly 
constrained, a final option for reducing the speed is chosen. This option is to initiate 
the abort maneuver sooner when the vehicle is traveling at speeds sufficiently slow to 
make the required turn. 

Figure 5.21 plots altitude as a function of velocity for the maximum attack angle 
a m ax = 90° case. In this figure various times of interest are denoted with the letters 
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A-I. At A the abort is initiated at 183 seconds after launch. From A to B, the vehicle 
climbs to as to reduce the downrange distance from the abort site. From B to C, 
the vehicle speed is reduced by flying with the maximum attack angle as shown in 
Figure 5.22. At C, the turn towards the launch site has been completed as shown 
by Figure 5.23. From C to D, the vehicle empties its tanks by increasing its throttle 
as illustrated by Figure 5.24. This results in the climb shown in Figure 5.21. At D, 
the tanks are empty and the glide to the landing strip begins. During B to C and E 
to F the normal force inequality constraint is active. The vehicle climbs and reduces 
its attack angle (see Figure 5.22) to prevent the violation of this constraint. From G 
to H the vehicle loses excess energy by traveling in the high drag transonic region. 
Finally, at I the landing conditions are satisfied. 



Figure 5.21: RTLS Max. Alpha = 90, Altitude vs Air Speed 
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Figure 5.22: RTLS Max. Alpha = 90, Attack Angle vs Time After Launch 



Figure 5.23: RTLS Max. Alpha = 90, Groundtrack 
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Figure 5.24: RTLS Max. Alpha = 90, Thrust Profile 


Figure 5.25 illustrates the heading time history. This figure shows that the heading 
angle is slowly modulated between —220° and —300° as the vehicle approaches the 
landing site. This maneuvering assists in the removal of excess energy as the vehicle 
flies through the high drag transonic region. 
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Figure 5.25: RTLS Max. Alpha = 90, Heading Angle History 

Figures 5.26-5.28 show the RTLS air speed, altitude and groundtrack profiles for 
several values of a max . For the 90° and 75° cases, the vehicle initiates the abort at 
a significantly higher speed than the 60° case. All three cases show that the vehicle 
initially climbs, but the 90° and 75° cases pull much larger attack angles to break off 
excess speed through increased drag, before beginning to turn. This results in the 
elongated groundtracks of Figure 5.28. In all three cases, the climb serves to keep the 
length of the glide back to the launch site to a manageable size, and to prevent the 
kinetic energy from becoming excessive. 
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Time After Launch, sec 


Figure 5.27: Maximum Abort Time RTLS Maneuvers : Altitude Profile 
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Figure 5.28: Maximum Abort Time RTLS Maneuvers : Ground-track 


Figure 5.29 illustrates the attack angle histories for the three cases examined. All 
three profiles show that the vehicle flies at maximum attack angles until the turn 
toward the landing site has been made. For the 90° and 75° cases, the vehicle to 
slows down by maximizing drag. In the 60° case, the vehicle burns off the excess 
propellant by climbing. 
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Figure 5.29: Maximum Abort Time RTLS Maneuvers : Attack Angle Profile 


The throttle histories for the same three maximum attack angles are shown in 
Figure 5.30. Note that in the 75 ° and 60° cases, the throttle fluctuates during the 
initial phase of the flight. This appears to be a fairly typical case of chattering 
control, associated with linearly appearing controls (thrust) during deceleration with 
a quadratic drag polar. This figure stresses the important trade-off that is made 
between slowing the vehicle for the turn and decreasing the mass of the vehicle. The 
vehicle is attempting to simultaneously command speed, mass flow and flight-path 
angle rates which can only be physically realized in a time averaging sense. 
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Figure 5.30: Maximum Abort Time RTLS Maneuvers : Throttle Profile 


For RTLS maneuvers initiated at specified abort times, the maximum and min- 
imum payload masses that could be returned to the launch site were determined. 
Single and two-engine-out scenarios were investigated. For each of the cases exam- 
ined, the maximum payload mass that could be returned to the launch site was 80.25 
klbm. This is the maximum payload mass that can be delivered to orbit by a vehicle 
with a full complement of engines. 
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Chapter 6 


Summary 


In this study, optimal control theory was employed to investigate the performance of 
abort to orbit (ATO) and return to launch site (RTLS) maneuvers for a seven engine, 
winged-body single-stage to orbit (SSTO) vehicle which lifts-off vertically and lands 
horizontally. The maneuvers were initiated when the vehicle suffered an engine failure, 
while ascending to orbit along a maximum mass to orbit trajectory. Single-engine 
and two-engine out scenarios were studied. Different measures of performance were 
used for ATO and RTLS maneuvers. The optimal control problems are solved in 
discretized form via a non-linear programming (NLP) method. 

For the ATO maneuvers, the mass delivered to orbit was selected as the perfor- 
mance measure. For the single-engine-out scenario, ATO maneuvers were successfully 
made from abort times ranging for 0 to 100 seconds after launch. The corresponding 
payload mass delivered to orbit ranged from 65 to 78 klbm. In general, the per- 
formance increased as the abort time increased. Exceptions to this generalization 
occured as the abort time was varied from 50 to 60 seconds after launch. During the 
time frame, the vehicle enters the high drag transonic mach region which limits the 
maneuverability of the vehicle. The loss of maneuverability restricts the ability of the 
vehicle to correct the ascent trajectory to compensate for the loss in thrust. 

Two-engine-out ATO maneuvers were performed for abort times ranging from 
50 to 100 seconds after launch. The corresponding payload mass delivered to orbit 
ranged from 64 to 74 klbm. For abort times prior to 50 seconds, the thrust to weight 
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ratio of the vehicle is insufficient to make the ascent to orbit. 

One of the performance measures used for the RTLS maneuvers was the time of 
the abort. By maximizing the abort time, the window of decision for determining 
which abort maneuver (ATO or RTLS) to perform is maximized. For a single-engine- 
out scenario, maximum abort times were found for maximum attack angles varying 
from 60° to 90°. The corresponding abort times ranged from 130 to 183 seconds after 
launch. The vehicle carried the nominal payload mass of 20 klbm. 

Because the vehicle has single-engine-out capability from the pad, the conclusion 
was made that RTLS maneuvers should only be made if the mass of the payload 
was larger than that which could be delivered to orbit with an ATO maneuver. This 
assumed that an RTLS maneuver could be made at the time of the engine failure. To 
check this assumption, RTLS maneuvers which maximized the payload to be returned 
to the launch site were studied. RTLS abort maneuvers were initiated at times ranging 
from 0 to 100 seconds for single-engine-out scenarios and from 50 to 100 seconds for 
the two-engine-out scenario. The studies showed that the maximum payload mass 
that can be delivered to orbit by a vehicle with a full complement of engines (80.25 
klbm) could also be returned to the launch site for each of the combinations of abort 
times and number of engine failures examined. 

To determine the entire range of payloads that could be returned to the launch 
site, RTLS maneuvers which minimized the payload mass to be returned to the launch 
site were, also studied. These studies showed that the minimum payload mass was 
zero for each of of the combinations of abort times and number of engine failures 
investigated. 

Additional work in the area of SSTO abort performance might include optimally 
shaping the ascent trajectory of a healthy vehicle such that the window of decision 
for determining which abort maneuver ATO or RTLS to perform is maximized for 
range of specified payload masses to be delivered to orbit. Secondly, the numerical 
approach used to produce these results might be improved by incorporating a more 
accurate integration scheme for determining the state equation constraints. 
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Appendix A 


Development of 3-D OF Equations 
of Motion for a SSTO Vehicle 


In this appendix, a Newtonian mechanics approach is used to develop the equations 
of motion for a rocket traveling over a rotating Earth. This development is highly- 
influenced by [10]. 

A.l Reference Frames 

The reference frames used in this development are discussed in this section. The 
origin, fundamental plane and fundamental direction of each frame is defined. Illus- 
trations showing how the different frames are related are included. 

A. 1.1 Earth-Centered Inertial Frame 

The origin of the inertial frame (I-frame) is at the center of the Earth. The funda- 
mental plane of this system is the equatorial plane. The fundamental direction is 
along the vector connecting the center of the Earth to the center of the Sun at vernal 
equinox. 
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Figure A.l: I-frame and R-frame 

A. 1.2 Earth-Centered Rotating Frame 

Like the I-frame, the origin of the rotating frame (R-frame) is at the center of the 
Earth. The I-frame and R-frame also share the same fundamental plane, the equato- 
rial plane. The fundamental direction of this plane is coincident with the Greenwich 
Meridian at all times in the fundamental plane. Figure A.l shows the relationship 
between the I-frame and the R-frame. The I x and R x axes are respectively the fun- 
damental directions of the I-frame and R-frame systems. The fundamental planes 
are I x — I y for the I-frame and R x — R y for the R-frame. The R-frame rotates with 
angular velocity, u> (rotation rate of the Earth). 

A. 1.3 Position Frame 

The origin of the position frame (P-frame) is at the center of mass of the SSTO 
vehicle. The fundamental plane is the local horizontal plane and the fundamental 
direction is due east. Figure A. 2 illustrates the location and orientation of the P- 
frame relative to the R-frame. The fundamental direction of the P-frame is the Pg 
axis. The fundamental plane is the P$ — P# plane. The vector, f, is the position 
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Figure A. 2: R-frame and P-frame 
vector of the center of mass of the SSTO vehicle. 

A. 1.4 Velocity Frame 

The origin of the velocity frame (S-frame) is at the center of mass of the SSTO vehicle. 
The normal of the fundamental plane is the velocity vector relative to the R-frame. 
The fundamental direction is in the direction of increasing heading angle. Figure 
A. 3 shows the relationship between the P-frame and the S-frame. The fundamental 
direction of the S-frame is the axis. The fundamental plane of this frame is the 
5^ — S~, plane. 


A. 1.5 Body Frame 

The origin of the body frame (B-frame) is also at the center of mass of the SSTO 
vehicle. The normal of the fundamental plane is the vector pointing from the center of 
mass to the nose of the vehicle. The fundamental direction is directed along the right 
wing. In figure A.4 the relationship between the S-frame and B-frame is illustrated. 
The fundamental direction of the B-frame is the B y axis. The fundamental plane of 
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Figure A. 3: P-frame and S-frame 

this frame is the B y — B z plane. The vector, T, is the thrust vector which is directed 
along the longitudinal axis of the body, B x . 


A.2 Development of Equations of Motion 

As shown by figure A.l the R- frame rotates with angular velocity relative to the 
I-frame of 


ur/ i = u R z 

The angular velocity of the P-frame relative to the R-frame is 


(A.l) 


wp/R = 9 R z — 4> Pe 

The R z unit vector expressed in terms of P-frame unit vectors is 


(A.2) 


R z = sin <f> Pr + cos <f> 
Using relation A. 3 equation A.2 becomes 


(A.3) 
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S Y 

Figure A.4: S-frame and B-frame 

cDp/n — 9 sin <f> Pr — <j> P e + 9 cos <f> P^ (A.4) 

The angular velocity of the P-frame relative to the I-frame is found using 

Up/ 1 = U)p/R + LJr/I 

This expression yields 

&p/l = (9 + co) sin <f> P R - 4> P e + (9 + uj) cos <j> P# (A.6) 

The velocity of the Single-Stage To Orbit (SSTO) vehicle relative to the R-frame 
expressed in terms of P-frame unit vectors is 

f P = r Pr + r9 cos <j) P e + r<^> P (A. 7) 
In the S-frame this velocity is expressed as 

r s = v Sv (A. 8) 


(A.5) 
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The rotation matrix that rotates P-frame vectors to S-frame vectors is given by 


Sv 


cos 7 cos 0 cos 7 sin 0 sin 7 


Pe 

S $ 

= 

— sin 0 — cos 0 0 


h 

5 7 


— sin 7 cos 0 — sin 7 sin 0 cos 7 


Pr 


(A- 9 ) 


Using the inverse of this rotation matrix equation, A. 8 is expressed in terms of P-frame 
unit vectors 


rp = v sin 7 Pp + v cos 7 cos 0 Pg + v cos 7 sin 0 P# 
Equating expressions A. 7 and A. 10 yields 

9 = v cos 7 cos ^/(r cos <f>) 

<t> = vcos7sin0/r 
r = v sin 7 

j 

The angular velocity of the S-frame relative to the P-frame is 


(A.10) 


(A.11) 


tvs/p = 0 sin 7 Sv — 7 50 + 0 cos 7 5 7 (A. 12) 

Similar to equation A . 5 the angular velocity of the S-frame relative to the I-frame is 
found using 


uis/i = u) S /p -1- ujp/j (A. 13) 

Making use of expressions A. 9 and A. 6, equation A. 13 becomes 

us/i = 0sin7 + u/r(cos7sin7cos0tan0) -f w(cos0cos7sin0 + sin0sin7) 5y 
—7 -f v cos 7 /r + v cos (f> cos 0 50 

+0 cos 7 + v/r( cos 2 7 cos 0 tan 0) -f u;(sin <j> cos 7 — cos 0 sin 7 sin 0 ) 5 7 

(A.14) 

The position of the S-frame origin is 

7s<> = r sin 7 5y + r cos 7 5 7 (A. 15) 
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The absolute velocity of the S-frame origin (i.e. SSTO center of mass) is found using 

rso = r v S v + + r 7 S 1 + Co s/I x f So (A. 16) 

where 

ry = Component of S-frame Origin Along Sv 
r ij, = Component of S-frame Origin Along 
r 7 = Component of S-frame Origin Along 5\ 

After some simplification, equation A. 16 becomes 

rso = {v + u>r cos <f> cos 0 cos 7 ) Sv 

— (u;r cos <^sin^>) 5^ (A.17) 

— (tor cos <f> cos xj; sin 7) 

Similarily, the absolute acceleration of the S-frame origin is found using 

rso - v'vSv + v^Sip -f u 7 S 7 -f Hi ) S/I x f So (A. 18) 

where 

v v = Component of S-frame Origin Velocity Along Sv 
v ^ = Component of S-frame Origin Velocity Along S^ 
u 7 = Component of S-frame Origin Velocity Along S 7 

After a few pages of algebra, expression A. 18 becomes 
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rs<> = [u — w 2 r cos ^>(sin 7 cos <f> — cos 7 sin (j> sin ij>)] Sv 

+ to 2 r sin ^>cos <^>cos xf>/(ycosj) — 2u;(tan7Cos <^sin ip — sin <j>) 

+v cos 7 cos ij) tan <f>/r]v cos 7 S# 

+ [7 — u 2 r cos (j){ cos 7 cos <j> + sin 7 sin <j> sin ip)/v 
— 2cj cos (j) cos xl> — v cos 7/r] v 

(A- 19 ) 

The forces acting on the SSTO include : Earth gravitational force, SSTO thrust, 
aerodynamic lift and aerodynamic drag. The resultant of these forces expressed in 
the S-frame is 


F = (T cos a — D — mg sini) Sy 

+(T sin a sin fi + L sin g) S^, (A. 20 ) 

+(T sin a cos g + L cos g — mg cos 7) 5 7 
From Newton’s second law 


mr So = F (A. 21 ) 

This yields 

v = {T cos a — D)/m — g sin 7 + u 2 r cos ^(sin7cos <p — cos 7 sin 
ifr = (T sin a sin /z + L sin g)/(mv cos 7) — u 2 r sin cf> cos d>cos xf)/{y cos 7) 

+2o;(tan7cos fisinifr — sin^) — u cos 7 cos ^ tan <f)/r > 

7 = (T sin a cos g + L cos g)/(mv) + u 2 r cos <^>(cos 7 cos <j> + sin 7 sin <f> sin ip)/v 
+ 2 u> cos(f>cos xj) + ( v/r — g/v ) COS7 

(A. 22 ) 
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Appendix B 

NLP Scheme Used To Solve 
Optimal Control Problems 


The NLP scheme used to solve the optimal control problems in this study employs an 
existing optimization algorithm named NPSOL [8]. Developed at Stanford University, 
NPSOL is a set of Fortran subroutines designed to minimize a smooth function subject 
to constraints. NPSOL requires the user to supply a main or driver program and 
subroutines that define the objective and constraint functions and (optionally) their 
gradients. 

A group of subroutines that compose the multi2 library are used to convert optimal 
control problems into a NLP problems and supply NPSOL with the required objective 
and constraint functions subroutines. No gradient information is supplied by these 
subroutines. Therefore, NPSOL must calculate all gradients numerically. 

The user of this scheme supplies six subroutines and a driver routine. The six 
subroutines interface with the multi2 library to yield the objective and constraint 
functions. Figure B.l provides an illustration of the scheme. In this figure, the 
arrows point to a subroutine that is called by the routine at the end of the arrow. 
The user supplied subroutines and driver program are capitalized and the multi2 
library subroutines are represented with asterisks. 

A brief description of each of the routines in this scheme is given below. 
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Figure B.l: NLP Scheme 
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User Supplied Routines 


MADS3.f : Driver Program 

CDS3.f : Evaluates Cost Function 

BCDS3.f : Evaluates Boundary Condition Constraints 

PLTDS3.f : Evaluates State Equations 

CCNDS3.f : Evaluates State/Control Inequality Constraints 

SCNDS3.f : Evaluates State-Only Inequality Constraints 

PCNDS3.f : Evaluates Free Parameter Inequality Constraints 

Multi2 Library Routines 

: Supplies NPSOL with Objective Function Value 
: Supplies NPSOL with Values of Constraints 
: Evaluates Constraints Along Trajectory 
: Evaluates State Equation Constraints Along Trajectory 
: Evaluates State/Control Inequality Constraints Along Trajectory 
: Evaluates State- Only Inequality Constraints Along Trajectory 
: Calculates Nodal Midpoints 

NPSOL Routines 

npoptn.f : Supplies NPSOL with Optional Input Parameters 
npsol.f : Solves the NLP problem 

As authored by Dr. Dan Moerder, this scheme originally had only single phase 
capability. The author of this study extended the capabilities of this code to handle 
multiple phase problems. 


cstds3.f 

ds3con.f 
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